Quantized adiabatic charge pumping and resonant transmission 
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Adiabatically pumped charge, carried by non-interacting electrons through a quantum dot in a 
turnstile geometry, is studied as function of the strength of the two modulating potentials (related 
to the conductances of the two point-contacts to the leads) and of the phase shift between them. It 
is shown that the magnitude and sign of the pumped charge are determined by the relative position 
and orientation of the closed contour traversed by the system in the parameter plane, and the 
transmission peaks (or resonances) in that plane. Integer values (in units of the electronic charge e) 
of the pumped charge (per modulation period) are achieved when a transmission peak falls inside the 
pumping contour. The integer value is given by the winding number of the pumping contour: double 
winding in the same direction gives a charge of 2, while winding around two opposite branches of 
the transmission peaks or winding in opposite directions can give a charge close to zero. 
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I. INTRODUCTION AND SUMMARY 

Quantized|-adiabatic charge transport was first invoked 
by Thoulessja who considered the dc current induced in 
an infinite one-dimensional gas of non-interacting elec- 
trons by a periodic potential, which is also a slowly- 
varying, periodic function of time. The robustness of the 
quantization in such systems, with respect to the influ- 
ence of disorder or many-body interactions, was further 
discussed in Refs. || and [| Since then, the focus of i 
terest in this phenomenon, termed "electronic pump' 
has been shifted to investigations of confined nanostruc- 
tures, e.g., quantum dots or carbon nanotubes, where 
the realization of the periodic time-dependent potential 
is achieved by modulating gate, voltages applied to the 
structure, in a "turnstile!' formjj a or by coupling to sur- 
face acoustic waves .E3~E£l The pumping of charge in such 
systems is achieved solely by the application of the time- 
dependent potential, and it exists even when the system 
is un-biased otherwise. Moreover, the charge transferred 
during a single cycle is independent of the modulation 
period. However, this charge is not necessarily quantized. 

The quantization of the pumped charge, that is, the 
possibility to transmit an integral number of electrons 
per cycle through an un-biased system, has been realized 
in devices which are weakly coupled to the leads, in which 
the (small) conductances of the quantum point contacts 
separating the dot from the leads are modulated in time. 
In such systems, the quantization has been attributed 
to the Coulomb blockade, which quantizes the number 
of electrons on the device. However, the modulation of 
the point contact conductances allows, in principle, for 
the possibility to cross-over from pinched-off conditions, 
where the Coulomb blockade is effective, to the almost 
open-dot conditions. In the latter, Coulomb-blockade ef- 
fects are expected to play a minor role as opposed to those 
of quantum interference of the electronic wave function 
over the entire structure. Nonetheless, the possible quan- 



tization of the charge pumped during a cycle through 
a large, almost open qua ntum dot, with vanishing level 
spacing, has been relatedlHMLa to Coulomb interactions 
within the dot. Pumping in an open nanostructure was 
realized experimentally for a quantum dot whose shape 
has been controlled by oscillating gate voltages.cl There, 
the amplitude of the pumped signal, which is indepen- 
dent of the modulation frequency, is found to increase 
with the driving force, though no quantization has been 
detected. 

While electronic correlations in fabricated nano- 
devices definitely play a role, it is still of fundamental 
interest to study electronic pumping of non-interacting 
electrons resulting from quantum interference effectSy, ta 
explore the circumstances under which it is optimalJl3'[L£l 
In this context, it is especially useful to investigate sim- 
ple, tractable models, where it is possible to relate the 
parameters characterizing the nanostructure, notably the 
conductance, with those that govern the magnitude of the 
pumped charge. Thus by considering a small, strongly 
pinched quantum dot, which supports resonant transmis- 
sion, it has been shownO that when the Fermi energy in 
the leads aligns with the energy of the quasi-bound state 
in the dot, the charge pumped through it during each 
period of the modulation is close to a single electronic 
charge. The correlation between resonant states and en- 
hanced pumping has also been found in Ref. [l8| Finally, 
it was pointed outtJ that the pumped charge in a model 
describing carbon nanotubes can change sign as the am- 
plitude of the modulating potential is varied. 

In this paper we explore in detail the relation be- 
tween resonant transmission and quantized charge pump- 
ing. The results we obtain can be summarized generically 
as follows. Consider a quantum dot, connected by two 
single-channel point contacts to its external leads. The 
individual conductances of these point contacts, and 
X r , are controlled by split-gate voltages which are mod- 
ulated periodically in time. Hence, during each cycle the 
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system follows a closed curve in the Xi — X,. parameter 
plane, which we call the "pumping contour" . As the var- 
ious parameters (e.g. the modulation amplitude P and 
phase shift 0, see below) are varied, this pumping contour 
changes its shape and location in the parameter plane, 
forming a Lissajous curve. (This can be also achieved 
by varying the gate voltage on the dot.) On the other 
hand, there are lines in that parameter plane, along which 
the transmission of the quantum dot is large. These are 
the "resonance lines" of the quantum dot. We find that 
the magnitude of the pumped charge and its sign are 
intimately related to the manner by which the pumping 
contour encircles parts of the resonance lines, and partic- 
ularly the peaks where the transmission is equal to unity. 
The charge is quantized when a significant part of the res- 
onance line is trapped within the pumping contour. For 
example, when the contour goes around the resonance 
twice, in the same direction, the charge attains the value 
of |2| in units of the electronic charge e. The sign is de- 
termined by the sense of the pumping contour. Thus, 
as function of the modulating amplitude, the pumped 
charge can vary as depicted, for example, in Fig. or in 
Fig. [7] bellow. 
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oscillate in time with frequency lo, such that the modu- 
lation amplitude is P and the phase shift between the Ji 
modulation and that of J r is 2<j). Then, the point con- 
tact conductances are given (in dimensionless units) by 
Xi = Jf and X r = J 2 . Explicitly, 

Ji = Jl + P cos(wi + </>), Xi = jf , 

J r = J L + P cos(tut - (/)), X r = J r 2 . (1) 

We solve for the charge pumped through— the quantum 
dot using the adiabatic approximation,™!^ that is, as- 
suming that the frequency lo is smaller than any charac- 
teristic energy scale of the electrons. Concomitantly, we 
determine the resonance lines in the Xi — X r plane. We 
then show that the pumping is quantized (see Fig. |l|) 
for values of P for which the pumping contour encloses a 
significant part of a resonance line. 

In our model, the couplings of the quantum dot to 
the leads are modulated in time and hence Ji and J r 
can attain both negative and positive values. This re- 
flects a modulation of the potential shaping the dot: the 
tight-binding parameters Ji and J r , which are derived 
as integrals over the site "atomic" wave functions and 
the oscillating potential, can have both signs. The ex- 
treme modulation arises when we set Jl = 0, so that 
the hopping matrix elements which couple the "dot" to 
the "leads" oscillate in the range {— P,P}. The dimen- 
sionless conductances Xi and X r are then modulated in 
the range {0,P 2 }. The corresponding pumping contour is 
then a simple closed elliptic curve. As the magnitude of 
the average hopping Jl increases, the shape of the pump- 
ing contour becomes more complex, turning into a Lis- 
sajous curve which folds on itself. These complex pump- 
ing curves yield a rather rich behavior of the pumped 
charge as function of the modulation amplitude P (see 
e.g. Figs. 1 and 7). 



FIG. 1. The pumped charge, Q, in 
units of e, as function of the modulating amplitude, P, for 
N = 4, ka = O.OOlvr, Jd = J, Jl = 2J, e = and <f> = OAn 
(See text). 

The resonance lines- in the parameter plane can be 
found experimentallyED When one of the point contacts 
is pinched off and the other is open, the conductance of 
the quantum dot, G, is dominated by the conductance 
of the pinched-off point contact, say Xi. Measuring the 
dependence of the quantum dot conductance on the gate 
voltage Ui controlling that point contact will yield the re- 
lation between Xi and Ui. In a similar way, one finds the 
relation between X r and U r . Using these relations and 
measuring G as function of Ui and U r , the resonance lines 
can be determined. 

To derive the above result we consider a simple model 
for a quantum dot, employing the tight- binding descrip- 
tion. In this model, the quantum dot is coupled to semi- 
infinite ID leads by matrix elements — Ji and — J r , which 



II. TRANSMISSION AND PUMPING THROUGH 
A QUANTUM DOT 

Our description of the quantum dot is a generalization 
of the model by Ng and Lee.EHl Imagine the quantum 
dot to be connected to the electronic reservoirs by two 
ID chains of sites, whose on-site energies are assumed to 
vanish, and whose nearest-neighbor transfer amplitudes 
are denoted by —J. The energy of an electron of wave 
vector k moving on such a chain is 

E k = -2 J cos ka, (2) 

where a is the lattice constant. From now on we measure 
energies in units of J. The confined nanostructure, that 
is, the quantum dot, is modeled by a 'bunch' of tight- 
binding sites connected among themselves. For simplic- 
ity, we take the latter in the form of a finite chain of N 
sites, each having the on-site energy eo, with a nearest- 
neighbor transfer amplitude — Jd- The effects of on-site 
interactions may be considered as part of eo, within a 
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Hartree approximation. This finite chain is attached to 
the left-hand-side lead at site 1, with the matrix element 
— Je, and to the right-hand-side lead at site N, with the 
matrix element — J r [see Eq. (Eh]. 

The calculation of the pumped charge requires the 
knowledge of the instantaneous (that is, when the time 
is frozen) scattering solutions of tlin problem, and the in- 

rixBE" 



stantaneous scattering matrix.l 



This scattering ma- 



trix also yields the instantaneous transmission coefficient 
of the quantum dot for any pair of the parameters Xi and 
X r , and hence the resonance lines in the Xg — X r plane. 
In the present case, the instantaneous scattering solu- 
tions can be straightforwardly obtained. Let us denote 
by xl the instantaneous scattering state at time t, which 
is excited by an incoming wave from the left reservoir, 
of energy E k , and similarly x r is the scattering solution 
excited by a wave coming from the right. We then have 



X\{x) = A a . t 



r t e 



on the left lead, 



Xe( x ) = A jt t e tkx , on the right lead, 
where x = na, and similarly, 



(3) 



Xl(x) = A , 



e~ lkx + r' t e ikx 



, on the right lead, 



Xr{x) = A 0jr t' t e lkx , on the left lead. 



(4) 



To normalize the incoming waves such that they will 
carry a unit flux, we put 



Ao,e = A . 



1 



2J sin feet 



(5) 



In Eqs. (||) and (|j), t t — t' t , r t and r' t are the instan- 
taneous transmission and reflection amplitudes, respec- 
tively. In our model, those are given by 



t t = t' t = -e~ lk(N - 1)a J e J r J D sin qaM k , 



n 



e ika X e X r sin(N -l)qa- J D X f sin Nqa ) M k 



-i2Nka 



(e ika XeXr sin(7V - l)qa - J D X r sin Nqa^j M k . (6) 



Here we have introduced the notations 
2i sin ka 



Mu = 



D k = Jg sin(7V + l)qa - J u e ika {X e + X r ) sin Nqa 



e l2ka X e X r sin(7V - l)qa, 



(7) 



and scaled all energies in units of J. The wavevector q 
describes the propagation of the wave on the quantum 
dot, such that 



E k — £q = — 2Jd cos qa. 



(8) 



We begin the analysis by determining the resonance 
lines of the transmission for our model. From the re- 
sult for the transmission amplitude t t [see Eqs. (0)], it 
is readily found that the transmission coefficient is given 
by 



T = \ttf 



Z 2 + ( J D sin ka sin Nqa(X e - X,.)) 2 



(2 Jd sin ka sin qa) 2 Xi X r 



(9) 



with 



E 

Z = j£ sm(N + l)qa + J D sin Nqa(X e + X r ) 
+ X e X r sm{N~ l)qa. (10) 



Clearly, one has T=l when Xi = X r and Z=0. For 
N > 1, these equations give two points on the diago- 
nal in the X^— X r plane. It takes some algebra to show 
that the second equation, Z= 0, corresponds to maxima 
of T when either X^ or X r is varied while the other pa- 
rameter is kept fixed. These local maxima occur on two 
"resonance lines" , shown in Fig. ||. This figure also shows 
a contour plot of the transmission in the Xg— X r plane, 
for the same parameters as in Fig. [l]. For the chosen 
set of parameters, the transmission is quite flat around 
the maximum; The resonance lines appear as a "ridge" 
on that plateau, with T decreasing slowly as one moves 
away from the maximum points on the diagonal along 
the resonance lines, but quickly as one moves away from 
these lines. 
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FIG. 2. A contour plot of the transmission in the X^— X r 
plane, for the same parameters as in Fig. [l]. The resonance 
lines are the colored curves in the figure. 
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Next we calculate the charge, Q, pumped through the 
quantum dot during a single period of the modulation. 
Here we employ an expression which gives Q in terms of 
matrix elements of the temporal derivative of the scat- 
tering ptitential, V , between the instantaneous scattering 
states.ESEj At zero temperature the charge per period 
reads 



Q = ^fdt (xiivixt) - {x\\v\x\) 



(11) 



This expression can be shown to reproduce the Brouwer0 
formula£3 

In our model, the temporal derivative of the potential 

is 



V(n,n') = - J e ( 8 n ,i5n',a + S n ',A,o , 

— J r ( <5 n ,iV<W,./V+l + 8n',Nfin,N+lj i (12) 

and hence the charge pumped from left to right is 

Q = -3?{^-^^(j,[x**(0) X *(a) -xt (O)xi(o)" 
+ j r [x t r *(Na) X t r ((N + l)a) 

-xf(Na)x\((N + l)a)])}. (13) 

This expression requires the knowledge of the scattering 
solutions on the two ends of the quantum dot. These are 
given by using Eqs. @ and 



x \(a) = jA^(j ka + r t e~ lka ), 
xl(Na) = ^„,Ae d *, 
Xtia) = ±A , r tte- tka , 

, t { AT ~\ 1 A ( —ikNa , / ikNa\ 
Xr( Na ) = Y ,r \ * ) 



together with the results in Eq. (g) 
the pumped charge, Q, is given by 



(14) 

It then follows that 



„ e f , Jr> sin Nqa sin ka 

Q = s.r — m — 



(X(X r — X r Xej Jr>Ek sin Nqa 
+ (xj± r - X 2 r ±e^j sin(N - l)qa 
+ (x r - 5Q) Jl sin(iV + l)qa . 



(15) 



The temporal integration in this expression, with the 
time dependence as given by Eq. (]]]), is now done an- 
alytically, separating the four poles of the denominator 
in cosuit. The results, for a selected set of parameters, 
are shown in Fig. [l]. 



We now relate the values of Q, for representative val- 
ues of the pumping amplitude P, to the pumping contour, 
which is the closed curve the system traverses during a 
single cycle in the X^— X r plane. Those curves are por- 
trayed in Figs. | (P=l), | (P=3), | (P=5), and | (P=10). 
In each of the these figures, the left plate shows the con- 
tours of equal transmission (in black) , the resonance lines 
(in red) and that section of the pumping contour which 
lies near the transmission maxima (in blue). The right 
plate exhibits the full pumping contour. (The numbers 
on the axes in the right plate indicate the scale of the 
pumping contour.) We begin with the situation at P=l. 
At this value, as can be seen from Fig. [l], the charge 
pumped is vanishingly small. At this value, the pumping 
contour indeed encloses only a small portion of the reso- 
nance line, far away from the transmission maxima (see 
Fig. |). 





FIG. 3. The pumping curve in the X|— X r plane, for P=l. 
The other parameters are the same as in Fig. |l|. The right 
plate shows the full pumping contour; The left one depicts 
the contours of equal transmission (in black), the resonance 
lines (in red), and the pumping contour (in blue). 

Next, we consider the case P=3, for which the pumped 
charge is close to —1 (in units of e). Examining Fig. |J, 
it is seen that in this case the pumping contour encloses 
only the upper resonance line, including the transmission 
peak there. For this value of P, the pumping contour be- 
gins to fold on itself (since and X r are definitely posi- 
tive), forming a Lissajous curve in that parameter plane. 
This tendency is of course enhanced as the pumping am- 
plitude P increases. The charge remains close to —1 as 
long as the pumping contour remains between the two 
branches of the resonance line, as in Fig. 0. 
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FIG. 4. The same as Fig. |, but for P=3. The other pa- 
rameters are the same as in Fig. |l|. 



Moving to the case P=5 (Fig. ||), at which the pumped 
charge is small and positive (see Fig. Q), we find that now 
the pumping contour encircles both transmission peaks 
on the resonance lines; The separate contributions of the 
two resonances to the temporal integration in Eq. (|l5| ) 
almost cancel one another, leading to a rather tiny value 
for the charge. 





FIG. 5. The same as Fig. §, but for P=5. The other pa- 
rameters are the same as in Fig. H. 



Finally, we consider the case P=10. The left plate in 
Fig. ^| is similar to the left one in Fig. ||, and indeed, 
the absolute value of Q is close to the one in both cases. 
However, the sense in which the pumping contour encir- 
cles the resonance is reversed in the two cases, as can be 
gathered by examining the right plates in Figs. |] and 
Hence, the pumped charge changes sign. 




FIG. 7. The pumped charge, Q, in units of e, as function 
of the modulating amplitude, P. The parameters are the same 
as in Fig. [j], except that cf> = 0.057T. 



The next five figures portray the pumping contours re- 
lation to the resonance lines for representative values of 
the amplitude P, for the parameters of Fig. [7. In Fig. ^ 
we have P=l. Then the pumping contour encloses just a 
small part of the upper resonance line, and also touches 
the peak on the lower resonance line. Indeed, Q has an 
intermediate value near —0.5, decreasing to zero as P 
moves away from 1 (see Fig. ^). 





FIG. 8. The pumping curve in the Xf — X,- plane, for P=l. 
The other parameters are the same as in Fig. |^. The right 
plate shows the full pumping curve; The left one depicts the 
contours of equal transmission (in black), the resonance lines 
(in red), and the pumping contour (in blue). 



FIG. 6. The same as Fig. |, but for P=10. The other 
parameters are the same as in Fig. |IJ 

We now examine the pumping when the phase shift <p 
of the modulation [see Eq. ([j])] is changed. Instead of 
taking <f> = OAir, which was used to plot the previous 
figures, we now choose (j> = 0.05-7T, keeping all other pa- 
rameters as before. Although one might have thought 
that this smaller value will reduce the magnitude of the 
pumped charge, we now find that |Q| attains the value 
of 2, see Fig. 0. 



Increasing the amplitude to the value P=3 reveals that 
the pumping contour encloses both peaks on the reso- 
nance lines, and therefore their separate contributions 
almost cancel one another, leading to a tiny value of Q. 
Also, the pumping curve begins to fold on itself, giving 
rise to the small "bubble" close to the origin (see Fig. ||) . 
Following the increase of that bubble as P in enhanced 
leads to the situation in which the bubble encloses the 
lower resonance, and then the charge attains the value 1 
(see Fig. [7]). As the bubble increases further, capturing 
the two resonance lines, as shown in Fig. [l0| Q again 
becomes very small. 
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FIG. 9. The same as Fig. |, but for P=3. The other pa- 
rameters are the same as in Fig. \A 



FIG. 12. The same as Fig. |, but for P=12. The other 
parameters are the same as in Fig. m 



III. CONCLUSIONS 




FIG. 10. The same as Fig. |, but for P=5. The other 
parameters are the same as in Fig. m 



As P grows on, we reach the interesting situation, de- 
picted in Fig. [n] for P=8, in which the bubble encir- 
cles twice the upper resonance line, leading to a pumped 
charge very close to 2e|. 




FIG. 11. The same as Fig. §, but for P=8. The other 
parameters are the same as in Fig. m 



Finally, Fig. [l2| shows the pumping contour for P=12. 
The contour is seen to shift away from the resonance 
peaks, yielding a vanishingly small value for Q. 



We have calculated the charge passing through a small 
quantum dot whose point contact conductances are mod- 
ulated periodically in time. The calculation is carried out 
for non-interacting electrons, under the assumption that 
the modulation frequency is much smaller than the elec- 
tronic relaxation rates. We have found the conditions 
for obtaining integral values for the pumped charge: The 
contour traversed by the system in the parameter plane 
spanned by the pumping parameters (in our case, the 
point contact conductances) should encircle a significant 
portion of a resonance line (along which the transmission 
of the quantum contact is optimal) in that plane. The 
magnitude and the sign of the pumped charge are deter- 
mined by that portion, and by the direction along which 
the resonance line is encompassed. 

The reason for this observation can be traced back to 
the expression for the pumped charge, Eqs. (|ll]) and 
( |l5| ) . The contribution to the temporal integration comes 
mainly from the poles in the denominator in Eq. (|T5|). 
These same poles are also responsible for the resonant 
states of the nanostructure. jbat j s, for the maxima in 
the transmission coefficient In this way, we have 
obtained a topological description for the phenomenon of 
adiabatic charge pumping. One can now imagine more 
complex scenarios: including higher harmonics of w in the 
time dependence of the point contact conductances can 
create more complex Lissajoux contours, which might en- 
circle portions of the resonance lines more times, yielding 
higher quantized values of the pumped charge. 

Finally, it should be mentioned that the results pre- 
sented above are obtained at zero temperature. At finite 
temperatures, the expression for Q should be integrated 
over the the electron-energy Ek , with the Fermi function 
derivative — 9//9i?fc.E3'E2l Hence, upon the increase of the 
temperature, the pumped charge would be smeared and 
suppressed. It would be very interesting to check the 
above predictions in more complicated models (for ex- 
ample, when there are several levels on each of the sites 
forming the quantum dot), allowing for a richer struc- 
ture of the transmission in the parameter space, and, of 
course, in real systems. 
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